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ABSTRACT 

Cosmological N-body simulations indicate that the spherically-averaged density pro- 
files of cold dark matter halos are accurately described by Einasto profiles, where the 
logarithmic slope is a power-law of adjustable exponent, 7 = dlnp/dhir oc r a . The 
pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable reg- 
ularity, and are well approximated by simple power laws, Q(r) = p/<r 3 oc r~ x . We 
show that this is expected from dynamical equilibrium considerations, since Jeans' 
equations predict that the pscudo-phasc-space density profiles of Einasto halos should 
resemble power laws over a wide range of radii. For the values of a typical of CDM 
halos, the inner Q profiles of equilibrium halos deviate significantly from a power law 
only very close to the center, and simulations of extremely high-resolution would be 
needed to detect such deviations unambiguously. We use an ensemble of halos drawn 
from the Millcnnium-II simulation to study which of these two alternatives describe 
best the mass profile of CDM halos. Our analysis indicates that at the resolution of the 
best available simulations, both Einasto and power-law PPSD profiles (with adjustable 
exponents a and Xj respectively) provide equally acceptable fits to the simulations. A 
full account of the structure of CDM halos requires understanding how the shape pa- 
rameters that characterize departures from self-similarity, like a or x, are determined 
by evolutionary history, environment or initial conditions. 
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1 INTRODUCTION 

The large dynamic range probed by current simulations of 
structure formation allows for robust measurements of the 
internal structure of large samples of dark matter halos span- 
ning a wide range of masses. The mass profile of CDM halos 
holds particular interest, mainly because of its direct connec- 
tion with key observational probes of halo structure, such as 
disk galaxy rotation curves, gravitational lensing measure- 
ments, and more recently because of the possibility of ob- 
serving the dark matter directly through its self-annihilation 
signal, or in laboratory detectors. 



E-mail: aludlow@astro.uni-bonn.de 



iNavarro et all |l99rj . Il997l , hereafter, NFW) argued 
that the simple 2-parameter formula 
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could be scaled to provide a good fit to the density profiles of 
simulated halos. The two physical scaling parameters are r s , 
the radius where the logarithmic slope of the density profile, 
7 = dlogp/d logr, equals —2 (the isothermal value) and the 
characteristic density, p c . As discussed by NFW, the param- 
eters, r s and p c , are not independent but rather reflect the 
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More recently, it has become clear that the mass pro- 
files of CDM halos are not strictly self-similar, as orig- 
inally suggested by NFW, but that a third parameter 
is actually required to describe their shape accurately 
llNavarro et al J 12004 IMerritt et all 120051. [20061; ICao et alj 



|200S| ; IHavashi fc White! |200S| ; INavarro et alJl2010D . These 
studies have shown convincingly that profiles where the 
logarithmic slope is a simple power-law of radius, 7 = 
— 2(r/r_2) a , provide excellent fits to simulated halo pro- 
files when the shape parameter a is allowed to vary. This 
implies a density profile of the form 
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which we call the "Einasto profile" for short jEinastoll 19651 ). 
The Einasto and NFW parameters are simply related by 
r-2 — r s and p-2 = pc/i. 

A second result that has received widespread atten- 
tion has been the fact that the pseudo-phase-space density 
(PPSD) profiles of simulated dark matter halos follow simple 
power-laws with radius: 
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This was originally reported by iTavlor fc Navarrcl (|200lf ) , 
and has been confirmed by a number of subsequent stud- 
ies, albeit with some debate over the actu al value of 
the best-fit power-law exponent (see, e.g., Rasia et alj 



2004lDehnen fc McLaughlin 20051; iFaltenbacher et alj|2007l : 



Vass et all 120091 ; iLudlow et all 1201(f) . Although originally 



shown to hold for the total velocity dispersion profile, eq. ([3]) 
also holds whe n a is replaced with oy, the radial velocity 
dispe rsion (e.g.. |Pehnen fc McLaug hlin 20051; INavarro et alj 
|2010| ). (We hereafter use Q when referring to the total 
pseudo-phase-space density profile, and Q r for its radial ana- 
log, p/o 3 r .) 

As discussed by ITavlor fc Navarrol (|200ll ) , the power- 
law nature of Q does not fully determine the halo density 
profile. Indeed, a wide range of different density profiles are 
consistent with this constraint, even for spherically sym- 
metric, isotropic systems in dynamical equilibrium. Well- 
behaved solutions, however, have only two possible asymp- 
totic inner behaviours, one where the central density di- 
verges like p oc r 2x_6 (the singular isothermal sphere is an 
example, for \ ~ 2), and another "critical" solution where 
the cusp asymptotic slope is much shallower; 7 — » — 2x/5 
as r -s- 0. ITavlor fc Navarrol (|200ll ) showed that the "crit- 
ical" solution for \ — 1.875 (the value pre dicted by the 
self-s imilar secondary infall solution studied bv lBertschingerl 
Il985l ) closely resembles the NFW profile over the radial 
range resolved by their simulations. 

The inner asymptotic behaviour predicted by eqs. ([2]) 
and (jSJ for the density profile is therefore different (Einasto 
profiles have a finite-density core rather than a cusp), imply- 
ing that simulations should be able to discriminate between 
the two. Indeed, the Q (r) profiles of E inasto halos are not 
in general power-laws l|Ma et alj|2009i ). which implies that 
halos cannot satisfy both eqs. ([2]) and (|3]) at once and that 
departures from either Einasto fits or power-law Q profiles 
should become apparent in simulations of adequate resolu- 
tion. 



Deviations from power laws in the oute r Q profiles of 
simula ted halos have already been reported bv lLudlow et alj 
(|2010t ). who argue that such departures are actually ex- 
pected in regions near the virial boundary separating the 
inner equilibrium region from the unrelaxed infalling enve- 
lope of the halo. Whether similar deviations are also present 
in the inner regions of halos has not yet been studied in de- 
tail. 

We address these issues here using a sample of equilib- 
rium CDM halos selected from the Mi llennium-II (MS-II) 
simulation (|Bovlan-Kolchin et alj|2009l ). We begin in Sec- 
tion [2] with a brief description of our simulations and anal- 
ysis techniques. Our main results are presented in Section 
[3] density profiles are presented in Section f3.1l and Q r pro- 
files in Section ^. 21 Since both Einasto and power-law PPSD 
profiles are incomplete dynamical models unless the velocity 
anisotropy profile is specified we analyze this in Section \3. 31 
In Section [3.41 and 1 3.51 we compare the Einasto and power- 
law PPSD models and directly asses their ability to accu- 
rately describe the mass profiles of simulated halos. We end 
with a brief summary of our main conclusions in Section [3] 



2 NUMERICAL SIMULATIONS 
2.1 Cosmological model 

Our analysis uses group and cluster halos identified in 
the MS-II, a ~ 10 10 -particle cosmological simulation of 
the evolution of dark matter in a 100 h^ 1 Mpc box. The 
run adopted a standard ACDM cosmogony with the same 
parameters as the Millennium simulation presented by 
ISpringel et all (|2005h : fi M = 0.25, Ov = 1 - SIm = 0.75, 
n s = 1, eg = 0.9, and a Hubble constant Hq = H(z — 
0) = 100 h km s" 1 Mpc' 1 = 73 km s _1 Mpc~\ but has 
better mass and spat ial resolution. Interested rea ders may 
find further details in iBovlan-Kolchin et alj (|2009l ). 



2.2 Halo selection 

We base our results on a sample of well-resolved equilibrium 
halos in MS-II. We focus on halos where N200, the number 
of particles within the virial radiufl exceeds 5 x 10 5 . This 
corresponds to group and cluster halos with virial masses 
above ~ 3.44 x 10 12 /i _1 M Q . 

In addition we require our halos to sa t isfy t he set of 
relaxation criteria introduced bv lNeto et alj (|2007I ) in order 
to minimize the impact of transient departures from equi- 
librium on halo profiles. These criteria include limits on (i) 
the fraction of a halo's virial mass in self-bound substruc- 
tures, /sub = M aub (r < r 2 oo)/M 2 oo < 0.07; ii) the offset 
between a halo's center of mass and its true center (de- 
fined by the particle with the minimum potential energy), 
doff = |rcM — r ccn \/r2oo < 0.05, and iii) the virial ratio 
of kinetic to potential energies, 2X/|$| < 1.3. All dynam- 
ical quantities have been computed in the halo rest frame. 



1 The virial radius, r200i defines the mass of a halo, M200, as 
that of a sphere, centered at the potential minimum, whose 
enclosed density is 200 times the critical density for closure, 
p clit = 3H 2 /8ttG. V200 = (GAf 2 ooA20o) 1/2 is the halo's virial 
velocity. 
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With these selection criteria our sample consists of 440 ha- 
los, more than half of the total number of halos (790) in the 
same mass range. 

2.3 Analysis 

Our analysis deals with the spherically-averaged density, 
p, and velocity dispersion, a, profiles of each of our 440 
dark matter halos. Each profile is built out of 25 spher- 
ical shells equally spaced in log 10 r spanning the range 
r CO nv ^ t ^ ^200 . Here r con v is the "convergence radius" 
defined bv lPower et al.l(|2003l ). where circular velocities con- 
verge to better than ~10%. Each spherical shell is centered 
at the particle with the minimum potential energy, which 
we identify with the halo center. 

For each radial bin we estimate the dark matter mass- 
density by dividing the total mass of the shell by its volume; 
the velocity dispersion is defined such that a 2 is two times 
the specific kinetic energy in the shell, measured in the halo 
rest frame. Analogously, we compute the radial velocity dis- 
persion, a 2 , in terms of the kinetic energy in radial motions 
within each shell. The velocity anisotropy parameter is de- 
fined by j3 = 1 — a 2 ^ n /2a 2 , where a 2 ^ = a 2 — a 2 . The p 
and o r values in each radial shell are also used to estimate 
the pseudo-phase-space density profile, Q r (r) = pjo 3 . We 
choose to focus our analysis on the radial Q r profiles (rather 
than the total Q = p/o~ 3 ) as this simplifies the analysis of 
Jeans' equations, in which the oy and j3(r) terms separate. 



3 RESULTS 

3.1 Density profiles 

Figure [T] shows the spherically-averaged density profiles for 
halos in our sample, with the residuals from best-fit Einasto 
profiles shown in the middle panels. All profiles have been 
plotted from r2oo down to the innermost resolved radius, 
r conv . Densities have been multiplied by r 2 in order to en- 
hance the dynamic range of the graph and to highlight halo- 
to-halo differences. Einasto fits are limited to the radial 
range r conv < r < 3 r_2, so that their parameters are not un- 
duly influenced by regions where crossi ng times are long, an d 
dynamical equilibrium may not hold (|Ludlow et alJuOlOh . 
Radii have been scaled to r_2, the radius at which the loga- 
rithmic slope is equal to the isothermal value, —2, and den- 
sities by p_2 = p(r-2)- 

We assess the quality of fits to the density profiles using 
the following figure-of-merit function, 

Nbin 

j=l 

All halos shown in Fig. [1] have i/>min < O.f. 

In order to highlight the differences in profile shape we 
have split our halo sample into three subsamples, each con- 
taining an equal number of halos: halos with a ^ 0.161 are 
shown in red; green curves show those with 0.161 < a ^ 
0.195, and blue curves are used for the rest. We adopt the 
same colour coding and halo samples in subsequent plots. 
For illustration, we also show three Einasto profiles with 
different values of the shape parameter, corresponding to 
the mean value of each halo subsample: a = 0.132 (dashed 



curve), a = 0.178 (solid black curve), and a — 0.230 (dot- 
dashed curve). 

The small residuals from best Einasto fits (typically less 
than 10% throughout the fitted region; see middle panels of 
Fig. [1} confirm that eq. ((2} provides an excellent descrip- 
tion of the density profiles of CDM halos down to the inner- 
most resolved radius. Further, Figure Q] also demonstrates 
that there is genuine variation in profile shape from halo to 
halo: profiles of different halos "curve" differently, yielding 
best-fit a parameter values that range between roughly 0.1 
and 0.25. The mass profiles of CDM halos are therefore not 
strictly self-similar, and a is a structural "shape" parameter 
genuinely neede d to describe accura tely the mass profile of 
CDM halos (seeH avarro et al.ll2O10l . for a full discussion). 



3.2 Pseudo-phase-space density profiles 

Fig. [2] shows the (radial) pseudo-phase-space density pro- 
files, Q r , of the halos in our sample, grouped as in Fig. [T] 
In order to take out the halo mass dependence, Q r profiles 
have been scaled radially by r_2, and vertically by P-2/V-2, 
where 1/-2 = \J G p-2 r 2 _ 2 is the characteristic velocity im- 
plied by the scale parameters p-2 and r_2. Our results con- 
firm the striking power-law nature of Q r profiles previously 
reported in the literature. For reference, we also show the 
power law, Q r oc r -1 - 875 (shown by a black dotte d line) , 
predicted by the self-similar solution o f lBertschingeil <|l985h . 

In order to emphasize the differences between halos, we 
show in the middle panels of Fig. [2] the residuals from the 
Bertschinger power law normalized to the mean value of 
P-2/V-2 f° r a-U halos in each sample. Although the residuals 
are in general small, the "curving" shape of their radial de- 
pendence provides a strong indication that the Q r profiles of 
simulated CDM halos deviate systematically from a simple 
power-law of fixed ex ponent y. 

As discussed by iLudlow et alj (|2010r ). the outer (r 2> 
r-2) upturn in the residuals is most likely associated with 
the transition from the inner, relaxed parts, to the unrelaxed 
outer parts, where infalling material has not yet had time to 
phase-mix with the main body of the halo. Su ch an upturn 
is pre sent also in the self-similar solution of I Bertschinger] 
(|l985h . and may be a general feature of the outer Q r profiles 
of CDM halos. Because of this, we have chosen the radial 
range r conv < r < 3 r-2 for all the fits we report here. 

Interestingly, the inner Q r profiles (r ^ 0.1 r-2) also de- 
viate from the r -1 ' 875 power law in a way that clearly de- 
pends on a. Note, for example, that at the innermost point 
the residuals of the lowest-a halos (leftmost middle panel in 
Fig. [2} are substantially larger than those of the largest-a 
halos (rightmost middle panel). We have explicitly verified 
that this a-dependence is not caused by grouping halos of 
different mass, nor by differences in numerical resolution. It 
is also unlikely to be due to anisotropies in the velocity dis- 
tribution, since, as we show below, all halos in our sample 
are nearly isotropic at radii this close to the center. We have 
also checked that this a dependence is not specific to our 
choice of Q r ', the total pseudo-phase-space density profiles, 
Q = p/o- :i , follow similar trends with a. 

These results imply that no single power-law can re- 
produce the radial dependence of the pseudo-phase-space 
density; if Q r does indeed follow a power-law of radius, then 




Figure 1. Spherically-averaged density profiles of all halos in our sample. All profiles arc plotted over the radial range T'conv ^ T *C ^"200 • 
Radii have been scaled by r_2; densities by p—2 = p(r~2)- Density estimates have been multiplied by r 2 in order to increase the dynamic 
range of the plot so as to highlight differences between halos. The full sample of halos has been divided into three subsamplcs according to 
the best-fit value of the Einasto parameter a: red curves (left panels) correspond to halos with a ^ 0.161; green curves (middle panels) to 
halos with 0.161 < a ^ 0.195 and blue curves (right panels) to halos with a > 0.195. There are equal numbers of halos in each subsample. 
To illustrate the role of the shape parameter, a, Einasto profiles with a = 0.132 (dashed), 0.178 (solid) and 0.230 (dot-dashed) are shown 
in the top panels. The middle panels are residuals from the best-fit Einasto profile with adjustable a; bottom panels are the residuals 
from the best-fit power-law Q r "critical" model (see Sec. 13.41 ) . In all cases fits are carried out over the radial < r < 3 r_ 2 . 



the exponent \ must vary from halo to halo. We show this 
explicitly in the bottom panels of Figure [51 where we plot 
the residuals from the best-fit power-law when the exponent 
X is allowed to vary. The small residuals over the fitted ra- 
dial range r conv < r < 3r_2) indicate that a power-law with 
adjustable \ provides a remarkably accurate description of 
the inner Q r profiles of CDM halos. 

Do Einasto profiles provide a better description of the 
spherically averaged structure of CDM halos than power-law 
Q r profiles, or vice versa? Because of different dimensional- 
ity, we cannot compare directly the goodness of fits to the 
p and Q r profiles shown in Figs. [1] and [2] One way to make 
progress is to compute the PPSD profiles corresponding to 
Einasto halos, or, alternatively, to compute the density pro- 
files of power-law PPSD models and compare them with the 
simulations. This may be accomplished by assuming dynam- 
ical equilibrium and solving Jeans' equations to link the p 
and Q r profiles, a procedure that, however, requires an as- 
sumption regarding the radial dependence of the velocity 
anisotropy, /3(r). We turn our attention to that issue next. 



3.3 Velocity anisotropy 

Velocity anisotropy profiles for all halos in our sample are 
shown in Figure [3] after rescaling all radii to r_2. Halos 
have been grouped according to the value of the best-fit a 



parameter, as in Figure [T] Solid lines with error bars show 
the median j3{r) profile of each group and the associated 
one-sigma scatter. The velocity anisotropy profiles exhibit 
a characteristic shape: they are isotropic near the center, 
become increasingly radial with increasing distance, before 
leveling off or becoming less anisotropic in the outskirts. 

FigureUshows the logarithmic slope- velocity anisotropy 
(7 vs j3) relation for the median profiles of each halo sub- 
sample using the same colour coding as previous plots. 
The dotted line s h ows t he linear /3(7) relation proposed by 
lHansen fc Moorel (|2006l ). which reproduces our simulation 
data well within the halo scale radius r_2 (i.e., for 7 > —2). 

The data in Fig. [4] also suggests that the P(-y) relation 
deviates from the Hansen & Moore fit in the outer regions 
(where j3 tends to a constant rather than the steadily in- 
creasing radial anisotropy predicted by Hansen & Moore). 
This dependence is captured well by the function 

/3(7) = ^(l + erf(ln[(A7) 2 ])), (5) 

as shown by the dot-dashed curves in Figure [3] The values 
of /9oo and A depend weakly but systematically on a, and 
are listed in Table [T] The radial /3 profiles corresponding 
to eq. are also shown with dot-dashed curves in each 
panel of Figure[3] assuming that y(r) is given by the average 
Einasto profile of each grouping. 
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Figure 2. Radial pseudo-phase-space density profiles, Q r = p/crf, for all halos in our sample. As in Fig. [T] halos are subdivided into 
three separate subsamples according to the value of their best-fit Einasto parameter a. Radii have been scaled by r_2, densities by p_2 
and velocities by V—i = r —2' '^ ne mean profiles and one-sigma scatter are shown as solid lines with error bars. The dotted curves 

in the top panels show a Bcrtschingcr r~ 1875 power-law to guide the eye. The middle panels plot the residuals from the Bertschingcr 
law normalized to the mean value of p_2/v^ 2 f° r an na l° s i n each sample. Lower panels show the residuals from the best-fit power-law, 
r~ x , with free-floating exponent, \- Fits are limited to the radial range r conv < r < 3r_2. The a-dependence of the Q r profiles is clearly 
evident from the systematic differences in the shape of the residual curves shown in the middle panels. 




Figure 3. Velocity anisotropy profiles, /3 = 1 — a^ an /2a^, for all halos in our sample. Profiles have been grouped according to the best-fit 
Einasto parameter, as in Figurc[T] Solid black lines with error bars show the median anisotropy profile and one-sigma dispersion for each 
halo subsample. The dot-dashed curves shown in each panel show the radial dependence of /3 (r) expected from eq. ((5J assuming the 
average Einasto density profile; dotted curves show /3(r) expected from lHansen fc Moorj ||200(J ) for the same p(r). 
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Hansen & Moore 
eq. (5) 




-2.0 -1.5 
7 = dlnp/dlnr 



Figure 4. Median density slope-velocity anisotropy relation for 
the three groups of halos s hown in previous figure s. The linear 
(3 — 7 relation proposed by lHansen fc Moore] ll2006f ) is shown as 
a dotted line; a relation of the form given in eq. 1(5} is shown 
separately for each sample using dot-dashed curves. There is 
some evidence that the mean /? — 7 relation departs from the 
lHansen &c Moorei feOOd) result in the outer regions of our halos. 
See text for further discussion. 



3.4 Density profiles from power-law PPSD 
profiles 

Once the radial dependence of the velocity anisotropy has 
been specified, equilibrium density profiles consistent with 
a power-law Q r (r) profile may be obtained fro m Jeans' 
equation. Following iDehnen fc McLaughlin! (|2005l ). we write 
Jeans' equation as 



3 

'5 K - 



2 x 



3 
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2x/3 



1/3 



(6) 



where y = p/p-2; x = r/r-2; ft-2 = 47rG p_2r 2 _ 2 / 'cr,-2 is a 
measure of the velocity dispersion in units of the "natural" 
velocity scale of the halo at r_2, and we have assumed that 
Q r (r) is a power law of exponent \- The logarithmic slope, 
dlnp/dlnr is denoted by 7, and primes indicate derivatives 
with respect to In a;. Once /3(r) is specified, the solution set 
of eq. (|6]) for given \ m fully characterized by k_2. 

It is instructive to consider the special case /3(r) = 0, for 
which three central asymptotic behaviours are possible: (i) 
a steep central cusp, 7 — > 2\ — 6; (ii) a central "hole" where 
y(0) = 0; or (iii) a "critical" solution with a shallow central 
cusp, 7 — > —2x/5- The latter may be thought of as the 
limiting case where the radius of the central "hole" solution 
goes to zero, and corresponds to a maxi mally-mixed state 
for g iven halo binding energy and mass (|Tavlor fc Navarrol 
l200lh . 

Given that steep central cusps such a s those in (i) 
are fi rmly ruled out by the data (see, e.g., iNavarro et al.l 
l201Cf ) , the "critical" solutions (iii) are the only viable density 



profiles consistent with the power-law PPSDs constraint. 
Although these considerations are strictly true only for 
isotropic systems, similar conclusions apply once /3(r) is 
specified. This is especially true considering that halos are 
very nearly isotropic close to the center (see Sec. 13 - 3 P - We 
shall hereafter adopt the "critical" solutions (evaluated nu- 
merically) as the density profiles corresponding to a power- 
law PPSD profile for a given value of X- 

The left panel of Figure [S] shows (in blue) the "critical" 
density profiles for three different values of x, and compares 
them to Einasto profiles. The values of a of the three Einasto 
profiles shown have been chosen to match as closely as pos- 
sible the profiles corresponding to the PPSD models. 

Clearly, for every value of \ it is possible to choose a 
value of a that matches the resulting density structure over 
a very wide range in radius. Deviations are only seen in the 
very inner regions, at less than 0.1% of the scale radius and 
therefore well outside the convergence region of any pub- 
lished halo simulation to date. For example, the black solid 
line in Fig. [5] shows the density profile of the billion-particle 
Aq-A-1 halo, whose convergence radius is roughly 0.01 r~i 
|Navarro et al.| [2010"). An Einasto profile with a = 0.17 and 
a "critical" solution for \ — 1-911 match the profile of this 
halo indistinguishably well. 



3.5 PPSD profiles from Einasto profiles 

Pseudo-phase-space density profiles corresponding to 
Einasto models can also be obtained by solving Jeans' equa- 
tion, which may be written as follows: 



ding 



dlnl/? „ a 
7 + —rr^ + 2/3. 



d In r d In r 

where = V c 2 /a, 2 . With /3(r) set by eq. ([5]), the right hand 
side of this equation is fully specified by r_2, P-2, and a. 
Solutions may therefore be found after choosing the value of 
( r (or its logarithmic derivative) at some fiducial radius, such 
as £_2 = Cr(r_2). The shape of the Q r p rofile is dictated by 
f_2, and it is not generally a power law l|Ma et al.ll2009h . 

Again, insight may be gained by considering the limiting 
behaviour of the isotropic solutions. Because Einasto profiles 
have finite central densities, oy must approach a constant at 
the center. The central velocity dispersion may be found 
numerically, and turns out to be quite insensitive to (-2 
(for given p_2, r_2, and a), provided that the £_2 is greater 
than about a tenth. (All of our halos are comfortably in that 
regime; indeed, the median value of (-2 for all of our halos 
is 1.274.) The finite value of the central velocity dispersion 
and its near invariance with £_2 also imply that the shape 
of the inner velocity dispersion profile is quite insensitive to 
the actual value of (_2- 

Therefore there is, in practice, a unique Q r profile that 
corresponds to an Einasto profile of given a, which we iden- 
tify with the single solution that is well behaved at all radii. 
This may be found by setting dln^/dlnr = a at r = 00; 
for a = 0.175, for example, this implies C,-2 = 1.265. 

The right panel of Figure [5] shows (in red) the PPSD 
profiles of Einasto halos, for three different values of a. For 
a = 0.1 and 0.17 the corresponding PPSD profiles are very 
well approximated by power laws over the whole plotted ra- 
dial range. Only for larger values of a, such as 0.3, are clear 
deviations from a power law noticeable. Even in this case, 
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Figure 5. Left: Einasto density profiles compared with "critical" profiles computed from Jeans' equation assuming power-law PPSD 
models. Three different values of a and of \ are shown, matched to highlight the similarity in the density profiles over a large radial 
range. Right: Radial PPSD profiles of Einasto halos compared with power-law models. As in the left panel, the values of a and x °f the 
three curves have been chosen so as to highlight the similarity of each pair of curves over a wide radial range. In both panels the thick 
solid curve in black shows the profile corresponding to the billion-particle Aq-A-1 halo, the highest-resolution CDM halo published to 
date. At the resolution of this halo, or lower, both Einasto and power-law PPSD models provide equally acceptable descriptions of the 
spherically-averaged structure of CDM halos. 



however, these are only evident in regions well inside 1% of 
the scale radius r_2, and therefore outside the converged re- 
gion of the highest-resolution simulation of a CDM halo pub- 
lished to date, the Aq-A-1 halo (shown by a solid black line). 
Deciding whether power-law PPSD models match CDM ha- 
los better than Einasto profiles, or vice versa, seems to re- 
quire simulations of even better resolution than Aq-A-1. 

3.6 Power-law PPSD vs Einasto density profile 
fits 

The results above suggest that both Einasto profiles and 
power-law Q r models provide equally good representations 
of the spherically averaged structure of simulated CDM ha- 
los, in spite of making very different predictions for their 
asymptotic inner structure. 

This conclusion may be verified quantitatively by fitting 
the "critical" solutions introduced in Sec. !3.4l to the density 
profiles of all halos in our sample and comparing them with 
Einasto fits. Residuals from the best fits obtained after vary- 
ing x are shown in the bottom panels of Figure [T] and are 
quite clearly indistinguishable from those obtained by fitting 
Einasto laws with adjustable a (shown in the middle panels 
of the same figure). 

Further quantitative evidence is provided in Fig. [6] 
where we plot the figure of merit, i/Wn, of the best Einasto 
fits compared with that obtained from the best-fitting crit- 
ical solution. Open and filled symbols in this figure are 
used to denote cases where, respectively, either isotropic or 



anisotropic (i.e., /3 given by eq. ©) critical solutions have 
been used in the fitting procedure. The close resemblance of 
the results obtained with either assumption implies that our 
conclusions are largely insensitive to our assumption about 
the radial dependence of the velocity anisotropy. 

Finally, our results imply a strong correlation between 
the best-fitting Einasto's a and "critical" x parameters for 
any given halo. We show this in Fig. [7] which suggests that 
the two parameters are essentially equivalent, and linked by 
a simple linear relation, x — 2.1 — 1.16 a. At the resolution 
of current simulations, it does not seem possible to decide 
whether Einasto or power-law PPSDs provide a better de- 
scription of the spherically-averaged structure of simulated 
CDM halos. 



4 SUMMARY AND CONCLUSIONS 

We use a sample of well-resolved equilibrium systems se- 
lected from the MS-II simulation in order to study the 
density and pseudo-phase-space density (PPSD) profiles of 
CDM halos. In particular, we explore the relation between 
Einasto profiles (often used to parameterize p(r)) and power- 
law PPSD profiles, which are often used to construct equi- 
librium models of CDM halos. 

We solve Jeans' equations to show that the PPSD pro- 
files of Einasto halos are close to power-laws, and that, con- 
versely, the density profiles of power-law PPSD models can 
be fitted by Einasto laws over a wide radial range. The two 
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Figure 6. Figure of merit, Vmini °f the best density profile fits to 
halos in our sample. Values on the y-axis are determined from the 
best-fitting "critical" power-law Q r model; those in the abscissa 
correspond to the best-fitting Einasto profiles. In all cases, fits 
are limited to the radial range r conv < r < 3 r_2 and have the 
same number of degrees of freedom: two scale parameters, p_2 
and r_2, and one shape parameter, \ or a - Halos in the shaded 
grey area have been omitted from the analysis. 



Figure 7. The relation between best-fit shape parameters, x an d 
a, derived from fits to the spherically averaged density profiles 
using cither Einasto profiles or "critical" profiles corresponding to 
power-law PPSD models. All fits are carried out over the radial 
range r conv < r < 3r_2, and assume the anisotropy profile given 
by eq. 10, with /3oo = 0.36 and A = 0.62. A dotted line shows 
the best-fit linear relations, \ = 2.1 — 1.16 a. 



Table 1. Parameters describing the velocity anisotropy profiles 
(eq. ©) 



descriptions differ, however, in their inner asymptotic be- 
haviour, although significant differences are only expected 
at radii well inside 1% of the scale radius, r_2, and are there- 
fore beyond the reach of current simulations. 

Our analysis of the MS-II halo sample shows that, 
over the resolved radial range, Einasto and power-law 
PPSD models provide an equally good description of the 
spherically-averaged structure of simulated CDM halos, pro- 
vided that the Einasto parameter a and the power-law ex- 
ponent x are allowed to float freely when fitting their radial 
structure. The strong correlation between best-fit a and \ 
parameters implies that they constitute, in practice, equiva- 
lent measures of the shape of the mass profile. These results 
confirm earlier suggestions that halo structure is not strictly 
self-similar, and that a "shape" parameter that varies from 
halo to halo is needed to characterize fully the structure of 
CDM halos. 

The spherically-averaged mass profiles of equilibrium 
CDM halos thus seem to be fully specified by three parame- 
ters: two physical scalings, r_2 and p_2, and a shape param- 
eter, q or x- Where CDM halos lie in this three-dimensional 
parameter space is likely linked to each system's evolution- 
ary history. Exploring the details of these relations remains 
a pending issue, but one that can be profitably addressed 
through large-scale numerical efforts such as the Millennium 
Simulation series. We expect to report progress in this area 
in the near future. 



(a) 


/3co 


A 


0.132 


0.253 


0.720 


0.178 


0.313 


0.656 


0.230 


0.361 


0.618 



ACKNOWLEDGEMENTS 

The Millennium-II Simulation was carried out as part of 
the programme of the Virgo Consortium on the Regatta 
and VIP supercomputers at the Computing Centre of the 
Max-Planck Society in Garching. JFN acknowledges sup- 
port from the Canadian Institute for Advanced Research. 
CSF acknowledges a Royal Society Wolfson Research Merit 
Award. This work was supported in part by an STFC rolling 
grant to the ICC. 



This paper has been typeset from a TgX/ DTffjX file prepared 
by the author. 



REFERENCES 

Avila-Reese V., Firmani C, Klypin A., Kravtsov A. V., 

1999, MNRAS, 310, 527 
Bertschinger E., 1985, ApJS, 58, 39 



The Density and Pseudo-Phase-Space Density Profiles of CDM halos 9 



Boylan-Kolchin M., Springel V., White S. D. M., Jenkins 

A. , Lemson G., 2009, MNRAS, 398, 1150 

Bullock J. S., Kolatt T. S., Sigad Y., Somerville R. S., 
Kravtsov A. V., Klypin A. A., Primack J. R., Dekel A., 
2001, MNRAS, 321, 559 
Dehnen W., McLaughlin D. E., 2005, MNRAS, 363, 1057 
Einasto J., 1965, Trudy Inst. Astroz. Alma-Ata, 51, 87 
Eke V. R., Navarro J. F., Steinmetz M., 2001, ApJ, 554, 
114 

Faltenbacher A., Hoffman Y., Gottlober S., Yepes G., 2007, 

MNRAS, 376, 1327 
Gao L., Navarro J. F., Cole S., Frenk C. S., White S. D. M., 

Springel V., Jenkins A., Neto A. F., 2008, MNRAS, 387, 

536 

Hansen S. H., Moore B., 2006, New Astronomy, 11, 333 
Hayashi E., White S. D. M., 2008, MNRAS, 388, 2 
Jing Y. P., 2000, ApJ, 535, 30 

Klypin A., Kravtsov A. V., Bullock J. S., Primack J. R., 

2001, ApJ, 554, 903 
Kravtsov A. V., Klypin A. A., Khokhlov A. M., 1997, ApJS, 

111, 73 

Ludlow A. D., Navarro J. F., Springel V., Vogelsberger M., 
Wang J., White S. D. M., Jenkins A., Frenk C. S., 2010, 
MNRAS, 406, 137 
Ma C, Chang P., Zhang J., 2009, ArXiv e-prints 
Merritt D., Graham A. W., Moore B., Diemand J., Terzic 

B. , 2006, AJ, 132, 2685 

Merritt D., Navarro J. F., Ludlow A., Jenkins A., 2005, 

ApJL, 624, L85 
Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 

563 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 
493 

Navarro J. F., Hayashi E., Power C, Jenkins A. R., Frenk 

C. S., White S. D. M., Springel V., Stadel J., Quinn T. R., 
2004, MNRAS, 349, 1039 

Navarro J. F., Ludlow A., Springel V., Wang J., Vogels- 
berger M., White S. D. M., Jenkins A., Frenk C. S., Helmi 
A., 2010, MNRAS, 402, 21 

Neto A. F., Gao L., Bett P., Cole S., Navarro J. F., Frenk 
C. S., White S. D. M., Springel V., Jenkins A., 2007, MN- 
RAS, 381, 1450 

Power C, Navarro J. F., Jenkins A., Frenk C. S., White 
S. D. M., Springel V., Stadel J., Quinn T., 2003, MNRAS, 
338, 14 

Rasia E., Tormen G., Moscardini L., 2004, MNRAS, 351, 
237 

Springel V., White S. D. M., Jenkins A., Frenk C. S., 
Yoshida N., Gao L., Navarro J., Thacker R., Croton D., 
Helly J., Peacock J. A., Cole S., Thomas P., Couchman 
H., Evrard A., Colberg J., Pearce F., 2005, Nature, 435, 
629 

Taylor J. E., Navarro J. F., 2001, ApJ, 563, 483 
Vass I. M., Valluri M., Kravtsov A. V., Kazantzidis S., 
2009, MNRAS, 395, 1225 



